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Solution of the Dirac equation predicts that when an electron with non-zero orbital angular 
momentum propagates in a cylindrically symmetric potential, its spin and orbital degrees of freedom 
interact, causing the electron's phase velocity to depend on whether its spin and orbital angular 
momenta vectors are oriented parallel or anti-parallel with respect to each other. This spin-orbit 
splitting of the electronic dispersion curves can result in a rotation of the electron's spatial state 
in a manner controlled by the electron's own spin z-component value. These effects persist at non- 
relativistic velocities. To clarify the physical origin of this effect, we compare solutions of the Dirac 
equation to perturbative predictions of the Schrodinger-Pauli equation with a spin-orbit term, using 
the standard Foldy-Wouthuysen Hamiltonian. This clearly shows that the origin of the effect is the 
familiar relativistic spin-orbit interaction. 

PACS numbers: 03.65.Pm, 03.65. Ge 
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I. INTRODUCTION 



The physical consequences of the spin-orbit interaction 
(SOI) for an electron in a spherically symmetric central 
potential are well-known: the corrections to the bound- 
state eigen-energies depend on the projection of the elec- 
tron's spin angular momentum (SAM) onto its orbital 
angular momentum (0AM), S ■ L [l|. This energy split- 
ting contributes to the famous fine structure of the en- 
ergy states in the hydrogen atom. One can calculate 
it using either the exact solution of the Dirac equation 
[3, |3|, which includes SOI implicitly, or by perturbation 
theory using the Pauli-Schrodinger equation, after ex- 
plicitly adding a spin-orbit term in the Hamiltonian [J] 
proportional to S • L. For an electron traveling within a 
cylindrically symmetric potential of infinite length, the 
energy states are continuous rather than discrete. How- 
ever there do exist transversely bound states, and one 
might still expect the SOI to alter the properties of these 
states in some way. Surprisingly, this simple and ana- 
lytically solvable problem does not seem to have been 
considered previously in the literature. 

In this paper we solve the problem of an electron trav- 
eling down a cylindrically symmetric step-potential that 
is translationally invariant in the z direction (see Fig. 1). 
We derive the wavefunctions and dispersion relations con- 
necting the electron's energy and momentum. We find 
in the cylindrical case that the energy corrections to the 
transversely bound states are proportional to the product 
ame, where a and mi are quantum numbers correspond- 
ing to the z-components of the electron's spin vector S 
and 0AM vector L, respectively. This stands in contrast 
to the case of a central potential, where spherical sym- 
metry dictates the dependence of the energy splittings 
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FIG. 1: Two concentric cylindrical surfaces with nearly equal 
radii a and a + Sa. The inner (outer) cylinder is positively 
(negatively) charged, thereby giving rise to an approximately 
constant electric field pointing radially outward between the 
cylinders, as expressed in equation ([ij. The electric field is 
zero elsewhere. 



upon quantum numbers j, £, and s only, where j, £, and 
s correspond to the electron's total angular momentum, 
0AM, and SAM, respectively. Similarly to the spheri- 
cal case however, the cylindrical SOI arises only in the 
presence of an an inhomogeneous potential; this interac- 
tion is absent for an electron in free space, even when 
considering axially localized beam-like states. 

The relationship between the electron's energy and lon- 
gitudinal propagation constant is given by the disper- 
sion curves for the distinct transverse states. We cal- 
culate the splitting of the dispersion curves induced by 
the SOI via two methods, paralleling the two standard 
approaches to the spherically symmetric case discussed 
above. First, we employ first-order perturbation theory 



on the Pauli-Schrodinger equation after explicitly adding 
to the Hamiltonian the appropriate spin-orbit term. In 
contrast to the former spherical case, we find that the 
added term is proportional to the product of the z- 
components of the spin and 0AM operators, SzLz- In 
the second approach we find nonperturbative solutions 
of the Dirac equation for the cylindrical geometry. The 
two results for the SOI splitting are found to agree in 
the appropriate limit, thus confirming the validity of the 
Hamiltonian used for the perturbative theory. 

The splitting of the dispersion curves has the following 
meaning, apparently found here for the first time: for a 
given electron energy, the phase velocity of the electron 
depends on whether the quantum number a has equal or 
opposite sign as the quantum number mi. That is, they 
depend on whether Sz points parallel or anti-parallel in 
relation to L^. This coupling of a and mi has an in- 
teresting consequence: it implies that there exist stable 
electronic states whose transverse spatial wavefunctions 
rotate as they propagate down the cylinder, with the di- 
rection of rotation depending on the sign of a (see Fig. 4) . 
One can therefore in principle exploit this interaction to 
achieve spin-controlled manipulation of the spatial elec- 
tron wavefunction. 

This spin-dependent rotational effect occurs in two dis- 
tinct contexts (see Fig. 3): when the electron wavefunc- 
tion is a superposition of degenerate energy eigenstates 
with the same value of a but opposite values of mi, the 
rotation occurs as a function of z. Complementarily, 
when the electron is a superposition of degenerate eigen- 
states of the ^-component of linear momentum, while still 
having the same a and opposite mi, the rotation occurs 
as a function of time. The possibility of this latter type 
of rotation for photons was predicted in Q. Both of 
these effects are the result of a varying relative phase 
between the propagating parallel and anti-parallel eigen- 
states, which in turn originates from the SOI-induced 
corrections to the dispersion mentioned above. Although 
these phenomena arise from relativistic dynamics, they 
persist even for nonrelativistic velocities. 

We are not aware of electron experiments to date that 
are sensitive to the predicted SOI effects in cylindrical 
geometry. Semiconductor waveguides used for studying 
ballistic transport of low-temperature electrons are typ- 
ically rectangular in cross section, so OAM is not con- 
served. Electrons in linear accelerator beams do not typi- 
cally have transverse coherence areas as large as the beam 
area, so coherent quantum effects would not be observed. 
In fact, the present calculation was motivated by consid- 
ering the analogous problem of a single photon traveling 
in a cylindrical optical fiber, where analogous effects have 
been predicted [g| , [3] . Although in this work we consider 
in detail only the simple case of a step-potential, we ex- 
pect the aforementioned SOI splitting effects to persist 
in any inhomogeneous cylindrical potential that is trans- 
lationally invariant in the z direction. However, if the 
requirement of translation invariance is dropped, we ex- 
pect the SOI to manifest itself in a more complicated 



way, in analogy with predictions of SOI for photons in a 
cylindrical Bragg cavity Q. In a future paper, we will 
elucidate the electron-photon SOI analogy in detail. 

The remainder of this work is organized as follows: in 
sectionini we derive the SOI Hamiltonian using a heuris- 
tic classical model of a charged particle with a magnetic 
moment propagating in a cylindrical waveguide. In sec- 
tion lllll we quantize this Hamiltonian and employ pertur- 
bation theory, thereby deriving the aforementioned en- 
ergy and propagation constant splitting, as well as the 
spin-controlled spatial rotation effect on the wavefunc- 
tions. After starting from the quasi-relativistic Foldy- 
Wouthuysen representation of the Dirac equation, we ar- 
rive in section [TVl at the same SOI Hamiltonian obtained 
in the heuristic model. We also give more explicit ex- 
pressions for the first order energy and propagation con- 
stant corrections and corresponding rotation rate. We 
provide in section|V]the most rigorous perspective on the 
SOI by obtaining relativistic wavefunctions directly from 
the Dirac equation with a step-potential, thereby show- 
ing the equivalence of our results in the Dirac, Foldy- 
Wouthuysen, and heuristic pictures in the appropriate 
limits. We conclude this work in section IVTl bv discussing 
the physical origin of the SOI for electrons, and briefly 
comparing it to the analogous case of a photon propagat- 
ing in a step-index optical fiber. 



II. SPIN-ORBIT HAMILTONIAN 

Consider a cylindrically symmetric potential which can 
be modeled by two concentric cylindrical surfaces with 
nearly equal radii a and a + Sa (see Fig. 1). The inner 
cylinder is uniformly positively charged (as observed in 
the laboratory frame) , and the outer cylinder is uniformly 
negatively charged, in such a way that overall the waveg- 
uide is neutral. The electric field is zero inside the inner 
cylinder and outside the outer cylinder, but is nonzero 
(and approximately constant) in the region between the 
cylinders, such that 



E = £o-e {p) P ~ SoO {p) p 

p 



(1) 



where Q (p) = 9 {p — a) — 9 {p ~ {a + 6a)) with 9 being 
the Heaviside step function and p the radial distance in 
cylindrical coordinates, and where p is the radial unit 
vector. The approximation on the right hand side of ([T|) 
is valid in the regime where da <^ a. The magnetic field 
is zero everywhere in the laboratory frame. 

We are interested in the case of an electron travel- 
ing down the cylinder with magnetic moment p and 
nonzero orbital angular momentum z-component (OAM) 
Lz = PPrfi with respect to the cylinder axis. We also 
assume that the electron is moving paraxially with re- 
spect to the cylinder axis such that |pt| <C |Pz|, where 
■pz = Pz^ and Pt = PpP + P4,<t' are the electron's longitu- 
dinal and transverse momenta in cylindrical coordinates. 



respectively. We will show that when such an electron 
is present in the region with nonzero electric field, the 
electronic motion gives rise to a spin-orbit interaction be- 
tween its magnetic moment ^-component fiz and 0AM 

The standard theory of SOI is summarized in Q- The 
magnetic field in the (primed) rest frame of the electron 
is 



B' 



V ^ 

-7— X E ; 
c 



X E 



(2) 



where v is the electron velocity in the laboratory frame, 
and the Lorentz factor 7 « 1 for sufhciently low v, which 
we will assume throughout this section. Also in ^, 
we have employed Gaussian units, following [9I|. The 
presence of the electron's magnetic moment fl in such 
a field gives rise to a magnetic dipole interaction en- 
ergy H' = —jl ■ B'. After accounting for the relativis- 
tic Thomas precession effect [13], which effectively con- 
tributes a factor of 1/2, this energy becomes 



H' 



i/i- 



X E 



2mc 



M 



(ex(p,z + pt)) (3) 



where p ^ PzZ + pt is the electron momentum in the 
laboratory frame. The SOI Hamiltonian therefore con- 
tains two parts in our present case with respective forms 
/I • (E X pzz) and /i • (E x pt). As the former term de- 
pends on the longitudinal momentum pz only, and there- 
fore does not involve the electron's transverse 0AM, we 
henceforth disregard it as a candidate for SOI. Upon 
employing ([1]), however, it is evident that the latter 
term involves a magnetic field vector proportional to 
Expr — Sqp^Q (p) z, which points either parallel or anti- 
parallel with the z-axis according to the sign of p^ = -Lz 
(see Fig. 2). From ([3]), this results in a SOI energy con- 
tribution of 



H: 



SOI 



1 £0 

2mc a 



lizLzQ (p) 



(4) 



where p k, a has been used. From ^ we see that when 
the electron is in the region a < p < a + 5a/it experiences 
a SOI energy shift proportional to the product of pz and 
Lz- In other words, the sign of the spin-orbit energy shift 
depends upon whether pz and Lz are pointing parallel or 
anti-parallel to each other. 



III. PROPAGATION CONSTANT SPLITTING 
AND SPIN-CONTROLLED ROTATION 

We quantize (g]) by letting pz -^ -i^Sz = -^^z 
and Lz ^ j-^ {oz is the Pauli matrix), so that the 
quantized Hamiltonian is 




B = Exp^ 



FIG. 2; The magnetic field contribution due to an electron 
propagating paraxially between the cylinders of the waveguide 
with nonzero p^, as experienced in the electron's rest frame. 
As discussed in the main text of the paper, we ignore the 
contribution due to p^ (represented by the dotted arrow in the 
figure), so that the field shown in the figure is that due only 
to the transverse component of momentum pr (represented 
by the bold arrow in the figure) . This effective magnetic field 
points in the negative- 2 direction for anti-clockwise p^ (as 
shown above), and in the positive-z direction for clockwise 
p^. 



H. 



SOI 




0(P) (5) 



where e = |e| is the elementary charge. 

The Hamiltonian in ([5]) is analogous to that which 
arises from an electron orbiting around a proton in a 
hydrogen atom — the canonical example for SOI. In that 
case, the electric field can be written as E = IgCouiombj.^ 
where gCouiomb ^ ^ is the Coulomb field due to the 
proton, so that the Hamiltonian in fS]) gives rise to the 
well-known atomic spin-orbit coupling Hamiltonian for a 
Coulomb potential: 



He. 



oulomb 



^0 



Coulomb 



2inn?c 



-S L 



(6) 



Though the Hamiltonians in ([5]) and ^ have similar 
forms and in both cases the SOI arises from the same 
Hamiltonian ([3]), the difference between the spherical 
and cylindrical geometries has significant physical con- 
sequences. In particular, for the cylinder case the spin 
and orbital quantum angular momentum operators cor- 
responding to the quantities pz and Lz commute with the 
Hamiltonian, while for the atomic interaction this is not 



the case, so that one must use the total angular momen- 
tum operator J^ and the z-component of total angular 
momentum J^ in the place of these. Therefore, while the 
total angular momentum quantum numbers j and rrij are 
good quantum numbers for the hydrogen atom, the spin 
and 0AM quantum numbers a and me are not. Con- 
versely, a, me and ttlj are all good quantum numbers for 
the cylinder case (though j is not, due to the breaking of 
the spherical symmetry) , so that states with well-defined 
a and me are energy eigenstates. We will make implicit 
use of this fact shortly. 

We treat (O as a perturbation of the standard 

Schrodinger Hamiltonian Hq = ^ — eV (p), where 
V [p) > 0. Our present task is therefore to find the 
unperturbed Schrodinger wavefunctions. Assuming the 

traveling wave form ip oc e'^^°^~~'i~ -' for the unperturbed 
eigenstates, in cylindrical coordinates the unperturbed 
equation of motion HqiJ) = Eotp takes the form 



SE„ 



ami 



ami 



Am?c'^ 



2ttN^- 



a-\-5a 



pdpjfmA (i^QP) 



2^2 ^2 ^ -^ImA i^oa) 



(10) 



Therefore, if the electron's SAM points parallel to its 
0AM mi, then the energy will shift upward, while for 
the anti-parallel case the shift will be downward. 

As introduced previously, two physical consequences 
of pO]) are the splitting of the phase velocity (and there- 
fore also the propagation constant /3o) of electron cylin- 
der wavefunctions with different values of cmie, and the 
related spin-controlled spatial rotation of these wave- 
functions. In order to better understand these effects, 
we note that due to the electron's wavelike properties, 
we can think of ([5]) as a dispersion relation defining 



v^v- 



KoV 



(7) 



where V^, is the transverse Laplacian V^ — -§t^ , and the 



transverse wavenumber is 



aP 



2 _ 






eV{p))-p'o 



(8) 



For a constant electric potential V [p) — Vb inside the 
cylinder, this is Bessel's equation, with solutions 



IV'o) - Af J|™,| («op) e^™'* (^ l'+_ ) e<P"''-^') (9) 



where we have constrained the wavefunctions to be fi- 
nite at the origin. In ^, N is a, normalization con- 
stant, the radial function J|,„j| (kqp) is a Bessel function 

of the first kind of order \mi\ = 0, 1,2, 



and 



5„- 



is a two component spinor composed of Kronecker delta 
functions such that 5a-\. = 1 if ct = -1-1 and 5^+ = if 
a — —\, etc. In expressing these wavefunctions, we have 
chosen the following complete set of commuting opera- 
tors, {H,pz, Lz, Sz}, which have the following respective 
eigenvalues, {Eq, H(3o, hme, fc}. We will henceforth des- 
ignate the states in ([9]) by |^o) = |»7i£,(T). 

From (O and © we conclude that the first-order cor- 
rection to the energy of an unperturbed state, SE = 
{mi, crlHsoilTni, a), is proportional to product ami pro- 
vided that the wavefunction is nonzero in the region 
a < p < a + 5a. This is indeed always the case for the 
transverse bound electronic states in ([9]) (we will show 
this in section [TVl when we apply the appropriate bound- 
ary conditions) . Explicitly, the first-order energy shift in 
this heuristic model is 



p{huj) 



2m 



{hhj + eVo) - k2 



(11) 



Later, in section |Vl we show that the Dirac boundary 
conditions imply in general that the value for an elec- 
tron's transverse wavenumber k differs slightly according 
to whether ami is positive or negative (that is, whether 
Sz and Lz are parallel or anti-parallel). We thus em- 
ploy positive and negative sign superscripts to denote 
these two cases, so that k, — > k"*" or k — > k^ depending 
on whether the the SAM and 0AM are parallel or anti- 
parallel, etc. Therefore, we conclude from pT|) that a 
parallel and anti-parallel state with the same value for 
P will have slightly differing frequency (energy) values 
Lo'^ and uj~ , respectively. This is the energy splitting 
which we have calculated in pH]) . However, we can also 
use PT|) to argue the converse — that parallel and anti- 
parallel states with the same frequency (energy) value uj 
will have slightly differing values for their propagation 
constants /3+ and j3~ , respectively. 

For a visualization of this point, refer to Fig. 3, which 
gives a plot of the dispersion relations (3 (E) for the states 
with ami = -|-1 and ami — —1, thereby explicitly show- 
ing the splitting of the curves (the dotted curve is a plot 
of the unperturbed dispersion relation). These parallel 
and anti-parallel states have different energies i?|| (E^) 
and E\y {E~) for a fixed value for the propagation con- 
stant /3o, as shown by the solid vertical and horizontal 
lines in the figure. Conversely, the two states have differ- 
ent (3 values /3|| (/3+) and /3[^ (/3~) for a fixed value of the 
energy Eq, as shown by the dotted vertical and horizon- 
tal lines. The horizontal and vertical arrows respectively 
show the directions (signs) of the energy and propaga- 
tion constant shifts 6E and SP for a parallel state (for 
an anti-parallel state, the signs of both SE and Sp are 
switched). The inlaid picture shows the resultant trans- 
verse spatial probability density distribution when the 




E„ 4.01 



E{MeV) 



FIG. 3: The splitting of the paraUel and anti-parallel states 
involving m = 1. This plot is a blown-up version of the inlaid 
box in Fig. 4, which plots the dispersion curves for all allowed 
states \mi, a) with R = 6 and eAV — 0.02mc^ (for definitions 
of R and AV, see equations (|27p and (|20p . respectively). The 
dashed curve is a plot of the unperturbed dispersion relation. 
The inlaid picture is a plot of the resulting transverse spatial 
probability density distribution when the parallel and anti- 
parallel states with \me\ — 1 are superposed. For a fixed /3, 
the azimuthal lobes of this distribution rotate as a function 
of time, while for fixed energy they rotate as a function of 
distance down the cylinder as shown in equation (|17p. In 
both cases the direction of rotation is dependent upon the 
spin of the superposition mode. For further discussion and 
the interpretation of the intersecting vertical and horizontal 
lines and arrows, see the main text of the paper. 



parallel (see equation p4|) ) and anti-parallel (see equa- 
tion P^ ) states with \me\ = 1 are superposed, as given 
by equation P7|) . 

In order to calculate the propagation constant shift 
6(3 to first order in terms of the energy shift SE — HSlu 
which we have already found, we expand the propagation 
constant f3 (huj) to first order in fko about the unper- 
turbed energy value ftwg, thus approximating /3 (hui) sa 
P{huJo) ±\SP\. In this way, the absolute value \S(3\ can 
be written as 



\Sf3. 



dp 



hduj 



h{uj — LOq) 



1 



hVz (uJo) 



\6E, 



(12) 



parallel state, SE^m^ is negative while SPame is positive. 
We therefore conclude that 



S(3a 






(13) 



This splitting in the propagation constants between 
electrons with parallel and anti-parallel spin and orbital 
angular momenta has a remarkable consequence: If one 
superposes a parallel and anti-parallel state with the 
same value for a and the same absolute value for mi, 
then the orbital angular momentum of the parallel state 
will be a jm^ |, while the orbital angular momentum of the 
anti-parallel state will be —<t \mi\. In the quasi-paraxial 
regime the resulting superposition wavefunction will then 
possess an azimuthal pattern that rotates as the particle 
propagates in the step-potential, with the sense of the 
rotation depending on the spin a. This spin-controlled 
rotation effect is a direct result of the varying relative 
phase between the parallel and anti-parallel states as they 
propagate down the cylindrical potential, which is in turn 
caused by the difference in the propagation constants of 
these states. 

More concretely, from (l9|), a parallel state has the ap- 
proximate form Q^ 

IV'ii) = ^J|™,| (^op) ( ll'l ) e<-l"^l^+/3^-)e-^^* (14) 
inside the cylinder, while an anti-parallel state is 



I'^lf) = ^-^Im,! [l^OP) 



5.- 



-i(dr|m£|0-^ z\-i-g-t 



where SE^m^ is the first-order energy shift, and v^ (wq) our final result, 

is the (positive) z-component of the group velocity of 

the matter wave, which is interpreted as the velocity of 

the electron as it travels down the cylinder. In order to 

determine the relative sign of 5(3ami and SEame, refer 

to Fig. 3. From the figure it is evident that for the 

dispersion curve of the parallel state, the energy shift 

6Ea-me IS positivc while the propagation constant shift 

SParuf is negative, as shown by the horizontal and vertical 

arrows. Conversely, for the dispersion curve of the anti- 



(15) 

The equal superposition of these two states, which we 
denote as IV'cr), is therefore equal to 

\^„) = ivj|„j,| {nop) ( £^ j 

(16) 

Note however that /e<'"l"^l'^+''^^) + e~*(''l"'l'^-'^"^)) 

can be written as cos {\mf\ -I- cr A/3z) e^^^ , where A/3 = 
\ {(3+ - /3~) and /3 = ^ (/?+ + /9") thereby leading us to 



^J|m,| i^op) cos (|m,| + aAPz) f f^+_ ) e<^^-^0 

(17) 

Equation (fT7|) is a major result of this paper; for 
a cylindrical step-potential it predicts the existence of 
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FIG. 4: Dispersion curves for each of the allowed transversely 
bound electronic states \mi,(j) for i? = 6 and eV = Q.02mc? . 
Note the splittings of the parallel and anti-parallel curves for 
states with \mt\ 7^ 0, which have been exaggerated by a fac- 
tor of 50 with respect to their actual values for purposes of 
visualization. The inlaid pictures are electronic transverse 
spatial probability density distributions associated with vari- 
ous values of \rai\, as discussed in the main text below. The 
varying azimuthal lobe structure of each of the plotted su- 
perposition states undergoes clockwise or counter-clockwise 
spin-controlled rotation as shown in equation (|17|) . A blown 
up plot of the small inlaid box which intersects the curves 
with \mii\ = 1 is presented as Fig. 3. 



Schrodinger wavefunctions with an azimuthal lobe struc- 
ture that rotate clockwise or counterclockwise about the 
cylinder axis as the particle propagates, with the sense of 
the rotation depending on the spin a. Furthermore, the 
rotation rate A/3 of the wavefunction has already been 
given implicitly via PU)) and P^ , and will be calculated 
explicitly using two different approaches in sections IIVI 
and El 

Dispersion curves associated with the parallel and anti- 
parallel states in equations (dH) and (fT5)) are plotted in 
Fig. 4 for several values of mp. The inlaid pictures as- 
sociated with each value of \mf,\ are plots of electronic 
transverse spatial probability density distributions. The 
distributions labeled by \mi\ = are obtained via equa- 
tion ([5]), while those labeled by \mi\ — 1, 2, and 3 are the 
rotating superposition states as given by (HZ]), which re- 
sult from superposing the parallel and anti-parallel states 
with equal \nrn\ from equations (fH)) and P^ . For a su- 
perposition state of a given energy, the spatial rotation 
rate A/3 is just half the splitting between its associated 
parallel and anti-parallel dispersion curves. The method 
used for calculating the curves is derived presently in sec- 
tion [IVl for a given value of \mi,\ we solve (^5)) and ([77]) 
for kq and then substitute the result into equation (|lip . 
For the values chosen for the figure, there are two allowed 
solutions for \m(,\ = 0, giving rise to two distinct \mi\ = 
dispersion curves and therefore also two distinct proba- 



bility densities. For \mi\ = 1,2, and 3, the dispersion 
curve splitting gives rise to stable superposition states as 
shown. 



IV. SPIN-ORBIT SHIFT: EXPLICIT 
CALCULATION 



In order to obtain explicit results for the SOI energy 
and propagation constant shifts, it is instructive to ap- 
proach the problem from the more rigorous viewpoint 
of the Foldy-Wouthuysen representation [J] , wherein the 
Dirac Hamiltonian has the general property that the pos- 
itive energy solutions are decoupled from the negative 
energy solutions so that we can describe the electron via 
a two-component spinor. In the presence of an arbitrary 
electrostatic field in the laboratory frame (the magnetic 

field is zero), to order (^) , the Dirac Hamiltonian in the 
Foldy-Wouthuysen representation takes the form Q 



Ha 



Zm 



8m^d 



■P' 



ieh 



2vn?c- 



■S • (E X p) 



■S • (V X E) 



V-E 



gj^2^2 



(18) 



where S is the spin vector operator of 2 x 2 Pauli matrices, 
the rest mass term has been dropped, and Gaussian units 
have again been employed. Our first goal is to argue that 
the contribution of the terms in curly brackets to the 
SOI has the form of the heuristically derived equation 
([5]). Note that the first term in curly brackets arises 
from the relativistic mass increase, and is independent 
of the form of the electric field E. In the canonical case 
of a Coulomb field, the next two terms (which are only 
Hermitian when taken together) give rise to the atomic 
spin-orbit interaction, while the last term becomes the 
well-known Darwin term. In light of this, we expect only 
the two middle terms to contribute to the SOI in the 
cylindrical case, and we henceforth drop the first and 
fourth terms. In section |V] we show that this is indeed 
justified by comparing the results of this section to those 
obtained directly from the Dirac equation. 

For electrostatic fields with zero curl the second term in 
the curly brackets also vanishes, so that after dropping 
the aforementioned terms there remains only the term 
2m^c^ S -(Exp). Furthermore, since E = —W, for a 
translationally invariant cylindrically symmetric poten- 
tial V {p) this term becomes 



2m 



2/2 



l dV{p) 
P dp 



S-(pxp) 



(19) 



(for a spherically symmetric atomic potential, this spin- 
orbit term has the same form, but with the replacement 
p — > r) . From p9p it is clear from the derivative that 



the spin-orbit interaction depends on the inhomogeneity 
of the potential and therefore does not occur in free space, 
thereby confirming the corresponding statements made in 
the Introduction. 

We now introduce the cyhndrical step-potential 



(20) 



V {p) = ^0 - ^ve [p - a) 



where Vb > and AV > 0. Substituting ^ into (fTO]) . 
and noting that — g'^ = S {p — a) (where S [p) denotes 
the Dirac delta function), we find that the Hamiltonian 
in (jISp takes the form (after dropping the aforementioned 
terms in curly brackets) 



where Hr, 



'0 - i;r 
Hamiltonian and 



H^=Ho + Hsoi (21) 

- eV (p) is the standard Schrodinger 



H 



AV 



SOI 



2m?c 



(5(p-a)S-(pxp) (22) 



is the perturbative SOI Hamiltonian. Furthermore, 
note that S ■ (p x p) in Hgoi can be expressed as 

SzLz + [ySx — xSyjPz- Since the unperturbed eigen- 

states \me, a) of Hq have already been given via ([5]) and 
([9]), we focus on the expectation of Hsoi as expressed 
in the unperturbed state basis, which is thereby propor- 
tional to the following two terms: 

{Hsoi) oi{m'f,a'\SzLz\mi>,a) 

+ {m'i,a'\[yS^-xSy\pz\mi,cF) (23) 

However, since {m'^^a'\{ySx ~ xSy\ Pz\'mt^^(7) always 

vanishes, we conclude that for the purposes of first-order 
perturbation theory we can write 



Hsoi — 



AV 



2rn?(? 



6{p- a) SzLz 



(24) 



Note that this is equivalent to our dropping of the term 
proportional to /i-(E x p^z) in ([3]). The operator Hsoi in 
([M)l is diagonal in the unperturbed basis, so we can read- 
ily calculate the energy shifts of the unperturbed eigen- 
states, 



5E„r 



TTh^eAV 



amg- 



N'JL, (Koa) 



(25) 



which agrees with the heuristically derived equation (|10p . 
since AV w EqSq. 

Though we have managed to obtain the general form 
of the energy shifts without considering the boundary 



conditions, we must do so now in order to obtain ex- 
plicit numerical results. We have already required that 
the wavefunctions be finite at the origin, resulting in ([9]), 
which is valid inside the cylinder. In addition to this, 
we furthermore constrain [V'o) to be zero at infinity, with 
both \ijJo) and its derivative continuous at the bound- 
ary (where p = a). For the region outside the cylin- 
der, the former condition results in the modification of 
^ via the replacement J\mi\ {i^op) -^ ^\me\ (i^op), where 
K\mt\ (^oP) is a modified Bessel function of the second 
kind of order \me\ (kq and kq denote the values of the 
transverse wavenumber inside and outside the boundary, 
respectively, as defined through ([8])). After emplo yin g 
the well-known cylinder function recursion relations [ll| , 
the latter two conditions thereby lead to the characteris- 
tic equation 



Ji 



Koa- 



mf 1 + 1 



(kqa) 



Ki 



J\niA ('^oa) 



KOO- 



raf 1+1 



(koo) 



i^i 



(Koa) 



(26) 



Equation ()26p is an equation in the two unknowns kq and 
Ko; in order to find a second equation in these variables, 
we use (PO)) to evaluate ([5]) inside and outside the cylinder 
and subtract the results to obtain 



«-S)«^-(T)X^)^«' 



(27) 



where A = h/mc is the electron's Compton wavelength 
[h is Planck's constant). Equations (|26p and (P7)) can be 
simultaneously solved for kq, and the result substituted 
into (O, which allows us to conclude that the wavefunc- 
tion is indeed nonzero at the boundary as required in 
section ini Finally, from ([9]), the normalization factor in 
^5t is found to be 



K\mt\-i (Koa)^|m,|+i ('«oa) 



N' 



™^ JLa ('*oa) 



J\me\~l (Koa) J\me\ + 1 («^oa) 



Jfm,\ ('^oa) 



(28) 



Therefore, (pS)) and P^ give the propagation constant 
corrections as 



SPam, = -anig 



1 heAV j J^|m<|-1 (^Qg) K\jne\ + 1 (^0«) 
J\me\-1 (^oo) J\mf\ + 1 (^oa) 1 



Vz 2m?c^a?' 



J? I (Kna) 



(29) 



so that A;3 = i ((S/3|| - 6(3^) = \ (^/9+|™,| - -5/3-1™,!) 
be written as 
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A/3 = 



\'mi\ 



1 helW \ K\„-,^\^i{Ki3a)K\^^\j^i{ki^a) 



Vz 2171^0^0? 



\mt,\ 



'\mi\ 



(koa) 
_i {Koa) J\mi\+i (koo) 



J? 



(Koa) 



(30) 



This is the exphcit form for the rotation rate of the elec- 
tron spatial wavefunction as defined in equation (J17p . 



such that its eigenvalue ctt = il, 



Lz + Sz is the total angular momentum operator and 
hr — 7573 Ip^T is the transverse helicity operator with 
1 
-1 
1 
-1 
while rrij is half an odd integer. 

For simplicity (and in order to avoid Klein's paradox 
as discussed below) , we will focus on the positive energy 
solutions to (l3T]l , which are of the form 



V. DIRAC EQUATION SOLUTIONS 

A few gaps persist so far in the development of this 
work. Specifically, in section IIIII we relied on the re- 
sult that the electron's transverse wavenumber k differs 
slightly according to whether ami is positive or nega- 
tive, while in section HVl we assumed that neither the rel- 
ativistic mass increase nor the Darwin term contributes 
to the SOI. Also, we have imphcitly assumed throughout 
the validity of the paraxial approximation, which is ex- 
pressed as |pt| ^ |Pz|, or equivalently as k <^ (3. In this 
section we will demonstrate the validity of each of these 
assumptions by deriving the relativistic analogue of equa- 
tions (|17p and ((30| . obtaining the bispinorial wavefunc- 
tions directly from the Dirac equation. Our derivation 
involves several steps. First, we construct the wavefunc- 
tions of interest, and boost them to a convenient frame. 
Next we apply appropriate boundary conditions and de- 
rive a characteristic equation. Finally, we approximate 
this equation to the appropriate order, thereby showing 
its equivalence to result ([50]) in section HVl 

The Dirac equation in bispinor form for an electron in 
a constant electric potential V {p) — Vq > is 



,2 
C(J ■ p 



mc ca ■ p \ I X 

2 



-mc 



X 



= i±E + eVo) 



(31) 



where E > is the absolute energy of the particle, and 
the upper and lower signs correspond to positive and neg- 
ative energy solutions, respectively. Free space solutions 
to the Dirac equation in cylindrical coordinates have 
been found [i2| ■ Since the potential V (p) is piecewise- 
constant for our case of interest, the solutions to ([31]) 
will have the same form (before boundary matching) as 



the ones in [I4I. Following [ij], we choose a complete 



set of commuting operators as {H,pT, Jz, hr}, with cor- 
responding eigenvalues {±E,hK,hmj,h<7T}, where Jz = 



\E, K, TOj, ±1) = 



Xnij 

hc{f3^iK) ^^ 

mc'-'+E+eVoXnij 



o<P-§t) 



(32) 



where 



1 (up) e 



^mj- 2 j0 



(33) 



±Z 1 (np) e 



and Zn (nop) denotes an arbitrary cylinder function of 
order n. Equation (j3ip contains the relativistic analogue 

of m. 



{cpf = h^c^ {f3^ + K^) ={E + eV {p)f - m^c^ (34) 

which for a sufficiently small step-potential [20|, can be 
used to derive a relativistic analogue to ([7f|) , 



(^-«>-^^(T) (S) 



vr^^u ,2 



(35) 



where E = ^m(? (in the laboratory frame) has been used, 

/ 2 N V2 . 

and 7 = I 1 — -^ 1 is the Lorentz transformation factor 
between the laboratory frame and the electron rest frame. 
Consider now the state 



Iff) = - (1^, K, mj, +1) + a\E, K, mj, -l)) (36) 

where mj = (m^ + icr) . By (|5^ and ([55)1 we can write 
this as 



W) 






Zmt (up) e 






mc-^+E+eVo 






=</3^-f0 



(37) 



r 



where 5a.± are Kronecker delta functions and a = ±1 
as before. We denote the state in ([57)1 by \a) because it 
becomes an eigenstate of S'z = -IS with eigenvalue fcr in 
the paraxial regime where k << (3. Note also that in the 
same limit \a) is also an eigenstate of L^ = —ih-^ with 
eigenvalue mi. 

It will simplify the analysis considerably to boost to a 
frame in which the terms involving (3 in both of the lower 
components of ([37]) become vanishingly small relative to 
those terms involving k. For an electron wave traveling 
with a sufficiently non-relativistic group velocity, such a 



frame will always exist provided that 



1 2 

■^mvq 



« eVo « 



mc 



where the lower bound ensures the existence 



of bound states, while the upper bound constrains the 
potential energy in order to avoid pair creation, which 
would invalidate the single particle Dirac theory. We 
henceforth assume that the above inequality holds and 
carry out the boost, so that in the new (barred) frame 
([57)1 is approximated as 



1^) 



0<7, + Zme 









(™^ + i)* 



2mc'^+7^eVo 



/ 



(38) 

where jz is the Lorentz transformation factor between 
the laboratory frame and the barred frame such that jz ~ 
7 since k << /3, and where Vq = jz^Vo « eVb and E « 
mc^ have been used. 

We now impose boundary conditions upon jo") by re- 
quiring the wavefunctions to be finite at the origin, 
zero at infinity, and continuous across the step-potential 
V (p) = Vo- Ave (p-a), similarly to section [IVl Note 
however that in the present case we drop the require- 
ment of the existence of a continuous derivative of the 
wavefunction at the boundary. The reason for this stems 
from the difference in order between the Schrodinger and 
Dirac cquations-the second order Schrodinger equation 
requires two conditions at the boundary (both continu- 
ity and a continuous derivative) in order to determine 
the wavefunction, while the first order Dirac equation re- 
quires only one. Application of these conditions on the 
boosted wavefunction (|38p results in the characteristic 
equation 



Jnii+a \HCl) 



2mc^ + 7zeVo Jme (na) 



Kme+a- (fia) 



2mc2 + jze {Vo ~ AV) K„r, (ka) 



(39) 



which, since 7zeVb << mc^, is well approximated by 



Jmt+cr (na) _ Kme+a (ko) 

K : K- 



Jjy 



Jmt \l^0') 



Kme (fia) 



eAV Jme+a (««) 

2mc^ Jmt (ko) 



(40) 

In Appendix A we show that equation (|40p is equivalent 
to the following condition: 



, J\nii\ + 1 {na) ^ ^|mf n-i («^a) 



J\mi\ (i^a) 



K\mt\ ('««) 



\mA 



zAV\ 



' 2mc^ 



"■ 






'-'Imfl {liC^) 



(41) 



Again, as in (j26p , we have arrived at an equation for two 
unknowns k and k, which is solved together with (j35p . 
The solution for k then yields /3 via (|M|) . 

From (|¥T|) we can clearly see that k (and therefore also 
(3) depends upon the quantity a-r^, which we assumed 
in section IIIII in order to arrive at the spin-dependent 



rotation effect of equation ^T7\ . In particular, cr-i^^ 



+1 in ([1T|) corresponds to the case of parallel spin and 
orbital angular nromenta (with k — > k+), while cr-r^ = 
— 1 corresponds to anti-parallel angular momenta (with 
K —> K~). In Appendix B, we show that equation (|4ip 
gives a prediction for the spatial wavefunction rotation 
rate 6(3 that agrees very well with that of equation (|30p . 
Therefore, we conclude that the Hermitian perturbation 



ieh 



S-{V xE)- 



2ni^c^ 



S -{Exp) (42) 



in (|T8l) is indeed the correct choice of Hamiltonian for the 
cylindrical spin-orbit interaction. For comparison, plots 

of the rotation rate A/3 as a function of i? = ^ J ^^Y 
as given by both the Dirac and Foldy-Wouthuysen ap- 
proaches are shown in Fig. 5. In both plots we keep 
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5>d0^° 


■ 


■ Dirac Solution 




4>d0'° 


■ 


A Foldy-Wouthuysen Solution 






■ 








■■■■■■■aaauHHB™, 











R 



FIG. 5: Plot of the Dirac and Foldy-Wouthuysen predictions 
for the rotation rate (propagation constant splitting) A/3 vs. 

for eV = 0.02mc^. 



R 



r 2Tva \ 



' eAV \ 



eV — 0.02toc^ constant, so that an increase in R corre- 
sponds to an increase in the ratio of the cyhnder potential 
radius to the Compton wavelength j. While it can be 
seen from the figure that the plots from both approaches 
agree very well with one another, there is a small rela- 
tive discrepancy which increases as R becomes small, as 
higher order relativistic effects come into play. Further- 
more, the predicted decrease in A/3 with increasing R is 
to be expected, as the transverse electron wavefunctions 
will tunnel into the step-potential with decreasing ampli- 
tude as J increases. The two predictions also approach 
each other asymptotically in this regime, as expected. 

Having demonstrated the equivalence of the Dirac and 
Foldy-Wouthuysen approaches with regard to the cylin- 
drical SOI phenomenon, our final aim is to derive the 
analogue of p?]) in section [1111 showing the spatial rota- 
tion of the Dirac bispinors. By an argument similar to 
that surrounding equation (J17p . starting from the (non- 
boosted) equation (|37[) we find that in the paraxial regime 
K « (3, for p < a, a parallel bispinor (that is, ami = +1) 
has the form 



l^p)=(±l)'"' 



-><y,+ 



mc^+E+eVo °'' + 

, -hc0+ c I 

\ mc^+E+eVo '^- / 



X J|„,,|(«;+p)e<^l'"^l*+''^^)e-'^' 



while an anti-parallel bispinor [am^ = —1) is 



(43) 



IV'ap)=(±l) 



\mi\ 



hc/3~ 
mc-'+E+eVo ' 



-;ic/3~ 



-S„- 



X J\ 



\ mc'^+E+eVo^ 

i[ij\mi\ip-p^ z) -i^t 



|mf| l« P) e 



(44) 



The key point here is that we are working in the near- 
paraxial regime, where k+ w k^ w S = 2 (.1^^ + k,~) 
and /3+ « /3- = /3 = i {13+ + P'). This fact allows us 
to make the following approximation: we completely ne- 
glect the small differences in transverse wavenumber n+ 
and propagation constant (3+ in the amplitudes of the 
spinorial components of (|43p and (j44[) . while retaining 
the propagation constant differences in the phase factors 
gi/3 z ^ This is a valid approximation, since a small vary- 
ing phase difference between propagating superposition 
states can have a large qualitative effect on the evolution 
of the probability distribution, while small amplitude dif- 
ferences will have only a small effect on this evolution. 
Under the aforementioned approximation, the approxi- 
mate superposition \^„) of (03]) and ((44|) can be written 
as 



IV;.) « (±1)1 



5 a.- 



J\mt\ (M 



gi((T|mf |0+/3+z) I -i(o-lmf|0-/3"z) 



(45) 



from 



section IIIII that 

can be written 



Recalling 

/i[(7\mi\4>-i-f3+z) _j_ ^-i{a-\me\ct>-(3^z)\ 

as COS {\mi\ (f) + a Af3z) e'^^'^, we present the final form 
for the quasi-paraxial spin-dependent spatially rotating 
Dirac bispinors: 



lV'a)=(±l) 



\mt\ 



\ 



mc2+_B-|-eVo '^■ + 



\ mc^ +_E+eVo 

X cos (|m^| -f aAl3z) e*('^ 



J\me\ (l^P) 



(46) 



If we neglect the two lower small components, we find 
that ([in]) does indeed approximately reduce to the two- 
component Schrodinger spinor in (|17p , and clearly shows 
the spatial rotation effect. 



VI. CONCLUSIONS 

We have shown via direct solution of the Dirac equa- 
tion for a cylindrical step-potential that the SOI Hamil- 
tonian derived heuristically as equation ([5]) and more 
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rigorously as equation p4|) correctly predicts a splitting 
of the dispersion curves of the electronic eigenstates ac- 
cording to the relative direction of their spin and orbital 
angular momenta. This splitting can cause a propaga- 
tion constant (phase velocity) difference between parallel 
and anti-parallel states, which in turn gives rise to sta- 
ble states that exhibit spin-controlled rotation of their 
spatial probability distributions. In particular, we found 
that for a given energy, a parallel electronic state has 
a slightly smaller propagation constant than that of an 
anti-parallel state. Although we have treated only the 
simple case of a step-potential in detail, it is clear from 
(|19p that any inhomogeneous cylindrical potential that is 
translationally invariant in the z direction will give rise 
to a similar spin-orbit interaction. 

Another way of looking at the difference between par- 
allel and anti-parallel states is the following: numerical 
solution of equation (f4T|) implies that a parallel elec- 
tronic state has a slightly larger value for its transverse 
wavenumber k as compared to an anti-parallel state. It 
follows from this that the transverse radial wavefunction 
associated with a parallel state does not penetrate as far 
into the step-potential as that associated with an anti- 
parallel state. 

A similar SOI effect occurring for a photon propagating 
paraxially in a cylindrically symmetric step- index optical 
fiber can also be viewed in the above manner. For the 
photonic case, the step-index in the dielectric medium 
plays the role of the step-potential, and the photon he- 
licity plays the role of the electron spin. Stable, spin- 
controlled, rotating photonic superposition states with 
field distributions similar to those shown in Fig. 4 oc- 
cur also for the photon case [7|, which arise from a sim- 
ilar splitting of the dispersion curves for parallel and 
anti-parallel photons. We note that the photonic spin- 
controlled rotational effect (called the optical Magnus ef- 
fect) was predicted in Q for a graded-index fiber with a 
parabolic profile and also for a step-index profile, how- 
ever to our knowledge analytic results for the step-index 
case have not been presented in the context of the wave 
theory of SOI for a photon. 

As is well known, the basis of the electronic SOI is the 
sum of two physical effects: the interaction of the elec- 
tron's magnetic moment with the magnetic field resulting 
from the electron's motion through an inhomogeneous 
potential, and the Thomas precession resulting from the 
electron's curvilinear path of travel due to this potential. 
It is interesting, however, that for the analogous case of 
a photon propagating in an inhomogeneous medium, the 
SOI effect persists although the photon lacks a physi- 
cal analogue to the electron's magnetic moment. The 
spin-orbit interaction of a particle with arbitrary spin has 
been discussed in [13], in which the SOI is explained in 
terms of non-commutative space-time coordinates which 
arise from a non-Abelian Berry gauge connection. In a 
future paper, we will give the details of the SOI calcu- 
lation for photons in a step-index fiber, employing both 
the "perturbative" and "exact" approaches in parallel 



with this present work, in order to further elucidate the 
electron-photon SOI analogy. 
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APPENDIX A: EQUIVALENT 
CHARACTERISTIC EQUATION 

It will be convenient in what follows to express the 
Kn (x) functions (modified Bessel functions of the sec- 
ond kind with purely real arguments) of (|40p in terms of 
Hn (ix) (Hankel functions of the first kind with purely 
imaginary arguments). These functions are related by 
the identity [Ij] 



Hn i^x) = (-Z)"+1 -X„ (X) 
TT 



(Al) 



The reason for this replacement is that the J„ (x) and 
Hn (ix) obey the same recursion relations J^] while the 
J„ (x) and Kn (x) do not. Using (|A1[) . we find that equa- 
tion (Uni) can be written as 



Jme+a (u) H,ne+a («) 

U V- 



Jrr 



J, 



'Jrui v^J 



Hme (v) 



eu- 



m^+a 



(u) 



^riii \'^) 



(A2) 



where u = na, v = ika, e = 7z 2mc^ ' ^^'^ where we have 
for later convenience multiplied both sides by a. 

Our next task is to explicitly account for the absolute 
sign of nil in equation (jA2[) . We therefore replace rrii -^ 
± \mi\ in (jA2[) . where the upper sign corresponds to case 
where mi > (positive 0AM) in the spinor in ([37|l . while 
the lower sign corresponds to m^ < (negative 0AM). 
Under this replacement, (jA2p becomes 



J±[ 



J±\me\ (u) 



(u) ^ g±|m,|+a (t") ^ ^_^ J±\me\+'7 ju) 



H. 



±|?7l£ 



[v) 



J±\nii\ (u) 



We can re-express this, however, using the Bessel func- 
tion relations [ly| Z_„ (x) = (— l)"Z„(x) (the Z- 
functions stand for either the J„ (x) or the Hn (x) cylin- 
der functions), so that (lASp becomes 



Ji 



7111 |dl(T 



(u) ^H\„i^\±^{v) ^ ^^J\rni\±a_{u^ 



Now, for (T = +1 (IA4I1 becomes 



£U- 



(A4) 



J\ 



|m^|±l 



{u) 



J\mt\ (u) H\mi\ {v) 



H\rr.\±A-) ^ J\^,.\±^i-) (A5) 



J\mi\ (u) 
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which consists of two distinct equations, corresponding 
to the choice of sign in the expression \mi\ ± 1 appearing 
in the cyHnder function arguments in the numerators. 
Focusing now on the case involving \mf\ — 1, where the 
left hand side of (IA5D is 



J\. 



(u) 



Hi 



?71^| — 1 



(v) 



J\mt\ (u) H\^^\ [v) 

we employ the following Bessel function identity 



(A6) 



J\ 



nii | — 1 



(u) 



H, 



\me\-l 



{v) 



Ju 



{u) 



H\ 



mi\ 



iv) 



\mt\+l 



(u) 



m 



\mi\ 



(v) 



J\mi\ {u) 



H, 



\mi\ 



(v) 



(A7) 



which can be proved by substituting the fundamental 
identities [13] ^n-i (x) = ^Z„ (x) - Z^+i (x) (again, 
Z stands for either J„ (x) or iJ„ (x)) into the left hand 
side of (|A7p . Using (|A7p in (|A6p allows us to write the 
two equations in (jASp as 



J\m(\ + 1 (u) H\mi 

u—h — ■ — r^ V 



(v) 



Ju 



{u) 



H\mt\ (v) 



±e /lr-l^l^"^ (A8) 



Ji 



\me 



(u) 



We note here that (|A8p is an expression of (|A4p in the 
case where ct = +1 (note also the new factor of ±1 on 
the right hand side of ([Mjl ). 

We now insert cr = — 1 into (|A4p . and carry out a 
simplification analogous to (|A5p - (|A8p . concluding that 
for cr = — 1, (jA4p is equivalent to 



Te/'r'^^/"^ (A9) 



J\niA + l{u) i?|mJ+i(u) 

u — r^ V- 



J\me\ (u) 



H\mi\ {v) 



J\me\ (u) 



Comparing (jASp and (jA9p , it is apparent that both equa- 
tions can be expressed simultaneously via 



J\rm\+i (u) Hi 
u — z ^-s V 



J\mi,\ (u) 



H 



,nA + l{-) ^^^^^■J]rnA±Au)_ ^^^^^ 



\mi 



(v) 



J\rn,\ (u) 



Recalling that the upper and lower signs in (jAlOp denote 
the cases of positive and negative rrii respectively, we 
replace the "±" notation on the right hand side with the 
equivalent expression t^^ . In this way, we arrive at our 
final form for the characteristic equation. 



J| 



|m«| + l 



{u) 



H 



7ne I + 1 



(v) 



Jl 



J\m,_\ (") 



Hu 



(v) 



me 

a-. — -:£U- 



mi\-\-a 



(u) 



\mi\ J|™,| (u) 

(All) 
which, in light of (lA), is just equation (|4T|) of section 

El 

APPENDIX B: WAVEFUNCTION ROTATION 
RATE 



Our aim here is to expand equation (|4ip in such a way 
that it gives rise to an equation of the form ([50)1 for the 
wavefunction rotation rate A/3, thus showing that the ex- 
act solution presented in section HVl is well approximated 
by the pcrturbativc approach of section IIIII We start 
therefore with equation (j4ip . written in the form given 
in appendix A as equation (jAlip . 



J| 



\mf\ + l 



(u) 



J\m,\ (u) 



H\me\+i (v) me 
v^-— ^--- = o- -eu- 



\nii\+(7 



(u) 



Hu 



(v) 



me 



J\m,\ {u) 



(Bl) 



recalling that u -^ u'^ 
u — > u~ = K~a and v - 



= n'^a and v 



depending on whether a-r^ = +1 or cr-r^ = —1, re- 
spectively. Since e = 72 #^^ is small, however, we 



have from (IBll) that m"^ « m" 



We can exploit this 
+1 



by adding/subtracting equation (jBip with <J-r^ 
to/from equation (|BT|1 with ct-t^ = -1. 
Adding the two cases of this equation gives 



(it+) + u 3 (u ) — w+H (u+) — V H (u ) 



Ima 1 + 1 



{U+) 



V|m,|-l(w ) 



J\nie\ (u+) 



J\-mt\ (U ) 



where SJ (a;) 



"d + 



,(^) 



H(a;) 






(B2) 



We now 



J| |(x) ' ^.■"^ — Hi 11 

multiply both sides of (|B2p by | and Taylor expand 
9 (u*) and J ^ ' ]± ■. about the point w = 5 (u+ + u^), 

to first order in Su, where (5w = ^ {u'^ ~w~) so that 
u^ — u i Su; we also Taylor expand K (v) about the 
point w = 5 (f^ + v^), to first order in 6v, where Sv = 
i {v^ — v^) so that v^ = V zt Sv. These substitutions 
result in the following equation: 



13 



m3 (u) - m (v) + USuf 3' (m) - {Svf H' {v)\ 



eu 



mi^\ — 1 



{U) 



J\mi\+li'^) - J\mt\~l{u) ^^{ j\m,\ + l{u)+J\ 



2J2 I (u) 



(B3) 



where the primes denote derivatives with respect to func- 
tional arguments. In (B3), the term in curly brackets on 
the left hand side is negligible because it is second order 
in (du)" and (Sv) , while the term in curly brackets on 
the right hand side is negligible because both e and Su 
are small quantities. Thus, to first order we have 



uQ (u) — wH (v) 



eu 



J\mi\ + 1 W — J\mi\-1 (w) 



'\mi\ 



(u) 



and since psp implies that 



(B4) 



(B5) 



we have in (|B4[) and (jB5[) two equations in the two un- 
known variables u and v. Therefore, upon substituting 
for V via (jBSp . (|B4p can be solved for u numerically. 
We now subtract the equations (|Bip . and find 



V = l^| R^ 



+ ^y+) 



u ^ 



u-^ {u-) - (w+H (?;+) - v-R {v-)) 



+ J\rnt\ + l{u^ 



J\nn\ {u+) 



J\m.,\ {u ) 



(B6) 



Multiplying (|B6|) by ^ as before, again Taylor expanding 
to first order in 5u and 5v about u and v, and neglecting 
quantities of order {5u) and e5u, we arrive at 



which can be used in (|B7|) to yield 



^w^S iu) + S' (u) - 4h (w) + H' (w) 

U V 



£ / J"|mf| + 1 (m) + J \nie\~l (") 

2 



•^Imfl (^f) 



(BIO) 



To simplify the term on the left hand side, we sub- 



stitute ^(x) = 'r'''^r/ and H (x) = Jr"'^^"'/ , while 



'I" 



again using (lA) and (5B) along with the cylinder func- 
tion relations [l^l Z'^ (x) — Zn-i (x) ~ f ^n+i (x) and 
Z;^i (x) = Z„ (x) - ^Zn+i (x) {Z stands for either 
Jn (x) or Hn (x)), in order to obtain 



is (u) -H 3' (u) - in (v) + H' (w) = 



r?i^ I ^ J 



72 



(S) 



(Bll) 



For the term on the right hand side of (jBlOp , we use ^^ 
^ J„ {x) = Jn+l {x) + Jn-l {x) SO that 



£ f J\mt\ + l ju) + J|mt|-1 (") 
2 V ^|m,| (S) 



Ime] 



(B12) 



(5m5 (w) -I- uS' (m) fe - (5uH (u) -t- wN' (w) (5u) 






Ji 



\mt 



(u) 



(B7) 



Now, using (|55l) to expand v to first order, we find that 

uSu 



Rl± 



V ± 6v 



Rl 



(B8) 



so that 



5v 



u^ - i?^ 



= — ^u 



(B9) 



Substituting the resuhs ((BTT|1 and (|BT2)) in ((BT0| and 
solving for uSu then gives 



uSu « e |m^| 



\mi\ \ ' 

J\mi\-1 (u) J\mi\ + 1 (u) 






(u) 



(B13) 



Su can thereby be found by substitution of u and v as 
given numerically by (jB4p and (jB5[) . 

Having calculated the difference between the trans- 
verse wavenumbers Su = ^a {k~^ — k~), we can now find 
the difference between the associated propagation con- 
stants, A/3 = i (/?+ — /3~ ) . To do this we start with 
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equation (j34p . which upon substitution of _E « ^^Tnc? is 
equivalent to 



Substituting (|B13p into (|B17p then yields the desired ex- 
pression for A/3, 



{P' 



C2fi2 



(7^ - 1) m^c'' + 2-i^mc^eV 



+ [ev {p)f ^ c^re {k^Y 



(B14) 



in the unbarred (laboratory) frame. Taking the square 
root of both sides and using (7^ — l) = (7z^) and 
7zeVo << mc^ in order to Taylor expand the radical 
then gives 



/3±«7. 



eV I ( c eV 



1 / hn^ 



2 \-fzmVz 



(B15) 



Therefore, we find that 



A/3 



1 h 



4 jzmvz 



(B16) 



However, since K^a = u zt Su, to first order in 6u (jB16[) 
is equivalent to 



A/3 



1 



Izfnvz a^ 



-u5u 



(B17) 



A/3 



m^ 



1 feAy j ^|m,|-l (-w)-ft^|m,| + l (-w) 



u, 2rn?c'^a? 



\me\ 



i-iv) 

J\mi\-1 (u) J\mi\ + 1 (u) 



JU (^) 



(B18) 



where e = 7z 2 ^ '^^^ been used. For clarity, we remind 
the reader that 



-IV = ^| R^ 



2 via (IB5D. where i?? ig 



defined in equation 
equation (jB4p . 



and that u is found by solving 



Comparing this first-order result of (JBISP to equation 
([30|) . we see that the two equations are of the same form. 
Furthermore, since £ << 1 and 7^ w 1, equation (|B4p is 
nearly equivalent to equation (j26p . so that u « Koa and 

We have therefore demonstrated that 



i?2 



Koa. 



the perturbative approach of section HVl is approximately 
equivalent to the direct approach of section |Vl 
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Strictly speaking, the wave numbers k^ and k^ of the 
parallel and anti-parallel states are not equal to each 
other or to kq as given in (|14p and (|15p . However, the 
approximation k"*" « k~ ~ k" is indeed justified as we 
show near the end of section |V] via solution of the Dirac 
equation. 

We use a small potential satisfying eVb << mc? in or- 
der to avoid the violation of particle number conserva- 
tion, called Klein's Paradox. See also the discussion above 
equation p8)). 



